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We examine devices constructed out of multilayered sandwiches of semi-infinite metal-barrier-semi-infinite 
metal, with the barrier tuned to lie near the quantum critical point of the Mott metal-insulator transition. By 
employing dynamical mean field theory, we are able to solve the many-body problem exactly (within the 
local approximation) and determine the density of states through the nanostructure and the charge trans- 
port perpendicular to the planes. We introduce a generalization of the Thouless energy that describes the 
crossover from tunneling to incoherent thermally activated transport. 
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1 Introduction Many new technological developments are anticipated over the coming years in the 
field of nanotechnology. There is a current interest in trying to incorporate strongly correlated materials 
into nanoscale devices, because strongly correlated systems often have interesting bulk properties that can 
be tuned by changing the pressure, temperature, chemical doping, etc. What is lesser known is how the 
strong electron correlations are modified when the bulk materials are confined on the quantum scale and 
attached to other (noncorrelated) materials, like normal metallic leads. In particular, we expect there to 
be a reorganization of the electronic states driven by a charge transfer, associated with the mismatch of 
the chemical potentials, and by the normal-state proximity effect of the metallic leads on the insulator, 
which will produce exponentially decaying states within the barrier of the Mott insulator, analogous to the 
superconducting proximity effect in normal metals. Here we adjust our system to be overall charge neutral, 
so we do not investigate the charge transfer; we only investigate the proximity effect. 

From a device standpoint, the simplest type of nanostructure to create is a multilayered structure, where 
stacks of planes of one material are topped by another material, and so on until a given device and het- 
erostructure is made. Recent advances in pulsed laser deposition and molecular beam epitaxy have allowed 
many complex structures to be grown, with interfaces between the different materials being well defined 
up to a few atomic layers. 

The question we will address in this contribution is how do the properties of semi-infinite metal-barrier- 
semi-infinite metal multilayered heterostructure vary with the barrier thickness when the barrier is a Mott 
insulator, tuned to lie just slightly above the Mott metal-insulator transition (i.e. very close to the quantum- 
critical point, but on the insulating side). We examine both single-particle properties like the density of 
states (DOS), including a proof that the local DOS in the barrier of a single-plane barrier has metallic 
behavior at strong coupling. We also investigate transport, and discuss a generalization of the Thouless 
energy that is appropriate for Mott-insulating systems and governs the crossover from tunneling to inco- 
herent transport. We employ dynamical mean field theory (DMFT) to perform these calculations, which 
relies on the local approximation for the self energy to be accurate in these three-dimensional inhomoge- 
neous systems. 
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Fig. 1 Schematic of the multilayered nanostructure, where we take a finite number of self-consistent metal planes 
(N m = 4 in the figure, but N m = 30 in our calculations), couple them on one end to a bulk semi-infinite metal, 
and on the other end to a barrier described by the Falicov-Kimball model (with Nb — 2 in the figure depicted by the 
cross-hatched planes; in our calculations Nb = N varies from 1 to 30). 

2 Formalism We describe the Mott insulator by the spinless Falicov-Kimball model 1 1 1 



where iy is a Hermitian hopping matrix, CZj is the Falicov-Kimball interaction, and Wi is a classical variable 
that equals one if there is a localized particle at site i and zero if there is no localized particle at site 
i (a chemical potential /j, is employed to adjust the conduction-electron concentration). Since we are 
considering multilayered heterostructures (see Fig.Q, we assume that the hopping matrix is translationally 
invariant within each plane, as well as the Falicov-Kimball interaction. For simplicity, we will take the 
lattice sites to lie on the sites of a simple cubic lattice, with tij — t for all nearest neighbors (we use t 
as our unit of energy), and we will take Ui = U for all lattice sites i that lie within the barrier plane. 
This choice assumes that the bare kinetic energy is the same for the metallic leads and for the barrier, and 
that the barrier is strongly correlated via the Falicov-Kimball interaction (renormalizing its bandstructure) 
within the barrier planes. We also work at half filling, with fj, = and (wi) = w\ = 1/2. In this case, 
the chemical potential has no temperature dependence, and the electronic charge remains homogeneous 
throughout the system. 

The dynamical mean field theory for inhomogeneous systems was originally worked out by Potthoff 
and Nolting 1 2 1 and developed for these particular heterostructures in another publication [ 3 1 . Here we just 
include the relevant summarizing formulas. The starting point is to note that the system has translational 
invariance in the two-dimensional planar direction, so we can Fourier transform from real space to mo- 
mentum space; all physical quantities we will be interested in here depend only on the two-dimensional 
bandstructure e 2d = — 2t[cosk !r + coskj,]. We let a Greek letter (a, (3, 7, ...) denote the z-component 
of each of the stacked planes. Then, because an electron with energy e 2d in the perpendicular direction, 
decouples from electrons with different perpendicular energy, the problem for the Green's function reduces 
to a quasi-one-dimensional problem, that can be solved with the renormalized perturbation expansion 1 4 1 . 




(1) 
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The result for the local retarded Green's function at plane a is 

GJoj) = ( de 2d p 2d (e) - . 9 . . „ . , . / ■< — — ^ (2) 

y ' J y ' L a (e 2d ,u;) + R a {e 2d ,Lu) - [w + |*-E a (w) - e 2d ] 

with p 2d the DOS of a two-dimensional tight-binding square lattice and E Q (w) the self energy at plane a. 
The functions R and L are determined via recursion relations: 

R a+n {e 2d ,u;) = u + fj, - E a+n (ui) - e 2d - \ 2d — ; (3) 



R a +n+l{(- 2d , to) ' 



L a ^ n (e 2d ,u) = uj + ii-Y Ia - n {uj)-e 2d - ' ; (4) 

L a - n -i(e 2d , uj) 

for n > 0. These recurrences are solved by starting at n — ±oo for the right or left recurrence, and then 
iterating in n. Of course, in real calculations, we must assume R a = Roo and L a = L_oo for all a up to a 
finite distance away from the interfaces with the barriers; we include 30 such self-consistent planes in the 
metallic leads (on each side of the barrier) in our calculations. We determine R^ (£_oo) by substituting 
Roc (L_oo) into both the left and right hand sides of Eq. [Eq. @], which produces a quadratic equation 
that is solved by 



1 



R oc (e 2d ^) = LU + fl ^ bu 2 lk{uj) 6 ±y{u J + p-nui k (oj)-e 2d ] 2 -4 = L- 00 (e 2d ,u J ), (5) 

with the sign of the square root chosen by analyticity or continuity. In both Eqs. and (|4), we see that 
whenever the imaginary part of R or L is positive, it remains positive in the recursion, implying stability; 
similarly, when the imaginary part is zero, we find the large root is stable, which is the physical root. Hence 
the recurrences are stable. 

Once we have determined the local Green's function on each plane, we can perform the DMFT calcula- 
tion to determine the local self energy on each plane [5 6 1. We start with Dyson's equation, which defines 
the effective medium for each plane 

G^ 1 (a;) = G- 1 (w) + E a ( W ). (6) 
The local Green's function for the ath plane satisfies 

G a (u)) = (1 — Wi) ; htCi 1 ; , (7) 

with w\ equal to the average filling of the localized particles [note that this above form is slightly different 
from the usual notation |6|, because we have made the theory particle-hole symmetric by the choice of 
the interaction in Eq. ([0, so that /.i = corresponds to half filling in the barrier region and in the ballistic 
metal leads]. Finally, the self energy is found from 

E a M = G - a V)-G-V). (8) 

The full DMFT algorithm begins by (i) making a choice for the self energy on each plane. Next, we (ii) 
use the left and right recurrences in Eqs. Q and @ along with the bulk values found in Eqs. Q and the 30 
self-consistently determined planes within the metal leads to calculate the local Green's function at each 
plane in the self-consistent region from Eq. Q. Once the local Green's function is known for each plane, 
we then (iii) extract the effective medium for each plane from Eq. l|6), (iv) determine the new local Green's 
function from Eq. 0, and (v) calculate the new self energy on each plane from Eq. (|8j- Then we iterate 
through steps (ii)-(v) until the calculations have converged. 
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3 Numerical results The first thing we consider is the local DOS on the barrier plane for a single-plane 
barrier (N = 1) in Fig. 0(a). Note how the DOS looks like a Mott insulator for large U (in the bulk, the 
Mott transition occurs at U xs 4.9), with an upper and lower Hubbard band forming, but there is substantial 
low-energy DOS coming from the interface localized states that are actually metallic in character (negative 
second derivative of the DOS with respect to uj). In addition, we show the Friedel-like oscillations that 
develop in the metallic leads as we move away from the barrier interface (for U = 5) in Fig. 0(b). Note 
how the oscillations have a shorter period, and a smaller amplitude as we move away from the interface. 
Already at 10 lattice spacings away from the interface, we see there is only a small difference from the 
bulk DOS (which is why 30 self-consistent planes is sufficient for this work). 




Frequency Frequency 



Fig. 2 (a) Local DOS for the barrier plane when TV = 1 and for various U (indicated by the numerical labels). 
The inset highlights the low-energy region, where one can see a low-weight metallic DOS form for large U from the 
interface localized states of the nanostructure. (b) The local DOS in four barrier planes for U = 5 and N — 5 (the 
first barrier plane is at a = 0; the metallic lead runs from a = — 1 to a = —30 on the left hand side). Note how the 
amplitude of the Friedel-like oscillations are quite small even by the time we hit the tenth plane from the barrier. 



We show that the peak value at w = is easy to derive from an analysis of the Potthoff-Nolting 
algorithm, and that it is nonvanishing for all U. To begin, we consider the non-self-consistent solution, 
where we set L a = and R a = Roo for all a except a = 0, where we have a Falicov-Kimball 

interaction. Then the local Green's function in the barrier is 

G Q=0 M = / dep 2d {e)- — 1 (9) 

where the sign of the square root is chosen by either analyticity, or continuity. If we assume Eo is large 
(which occurs in the Mott insulator) then we can expand Eq. (|9jl in a power series in inverse powers of E 
to give 

G a=0 (u)n-—— v +... (10) 

with s(u) = J p 2d {e)^j(uj + \U -e) 2 -4. Taking the value for G a=o from Eq. ( II PL and plugging it 
into the self-consistent DMFT algorithm described above, allows us to solve for the self energy directly, 
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with the result 

If/2 

SoM = 9 T77T ? x C 11 ) 

ZUJ + U ± S(U!) 

which is large for U ^> 1, consistent with our ansatz. If we perform the integral in the definition of s(u>), 
we find it satisfies s(uj) = (ui + iU) ■ 0.653 + i ■ 1.05. Substituting this result into the self energy, and then 
into the Green's function, and evaluating the DOS, finally yields 

r \ 1-05 

Po{uj) * .[(i.o^ + ic/T (12) 

we compare this result to the exact calculated result in Tabled One can see the agreement is excellent. 

Table 1 Comparison of the results from Eq. i 1 2t and the exact numerical DOS for a single-plane barrier 
and various U. 



u 


approximate result [Eq. (11211 


exact numerical result 


6 


0.033 


0.037 


8 


0.0195 


0.021 


12 


0.009 


0.009 


16 


0.005 


0.005 


20 


0.003 


0.003 



In Fig. |3] we plot a false color (grayscale) plot of the DOS of the nanostructure in the near critical 
region U — 5 for a moderately thick (N — 20) barrier. Note how the Friedel oscillations are most apparent 
in the center of the band in the metallic leads (upper part of the graph). We only plot the 40 planes on 
the left-hand-side of the nanostructure, because the symmetry of the structure guarantees the rest of the 
nanostructure can be determined by a mirror plane reflection. Note how there are few oscillations within 
the barrier itself, and how the DOS rapidly becomes small at low energy (because it is an insulator). It is 
possible to even see some oscillations induced in the metallic lead at energies close to the band edge, and 
at positions close to the interface. 

We examine transport properties in Fig. |4] The resistance of a junction is calculated in the linear- 
response regime via a Kubo-based formalism 1 7 1 with the current-current correlation function 1 8 3 1 . The 
formalism requires us to employ a conductivity matrix in real space, with matrix components correspond- 
ing to the z-axis label of the different planes in the system. The resistance-area product can be calculated 
for any temperature. 

We extract an energy scale from the resistance, which we call the Thouless energy 191 1101 . since it 
reduces to the well-known diffusive and ballistic limits, but it also defines an energy scale when the barrier 
is a Mott insulator. In the insulating phase, it is a function of time, and the point where Exh ~ T defines an 
important energy scale for the dynamics of the transport in a nanostructure — it signifies when the transport 
crosses over from tunneling behavior at low T to incoherent "Ohmic" transport at high T [see Fig. 0(a)]. 
On the other hand, the Thouless energy is inversely proportional to L, the thickness of the barrier for 
ballistic transport and inversely proportional to L 2 for diffusive transport [this can be seen from taking 
the low-temperature limit of Eq. Jl 31 and noting that the DOS is nearly constant near the Fermi energy at 
low temperature and recalling the R n oc L for diffusive transport and R n is independent of L for ballistic 
transport] . 

The generalized formula for the Thouless energy is II ll l3l 
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Fig. 3 False color (grayscale) plot of the local DOS for the near critical U = 5 nanostructure with N — 20 planes 
in the barrier. Note how the Friedel oscillations in the metallic lead are most apparent near the center of the band, and 
how there are limited oscillations in the barrier (the DOS decays exponentially fast with position in the barrier at low 
energy). (Color on-line.) 




0.01 0.1 1 0.01 0.1 1 

Temperature Temperature 

Fig. 4 (a) Thouless energy for the U = 5 nanostructure as a function of temperature on a log-log plot. The different 
curves are for different thicknesses of the barrier. The dashed line is the curve Em = T, and the special points 
that denote the crossover from tunneling to incoherent transport correspond to the points of intersection of the solid 
lines with the dashed line, (b) Resistance-area product plotted versus the temperature on a log-log plot. The different 
curves correspond to different barrier thicknesses. The solid dots and the dotted line plots the points where Eth = 
T . Note how the curves are flat at low temperature and for thin junctions indicating tunneling (but the tunneling 
resistance does not grow exponentially with the thickness when we are so close to the metal-insulator transition and 
at finite temperature-note the unequal spacing of the lines with N). At higher temperature, the resistance picks 
up strong temperature dependence, and the transport is best described by incoherent thermally activated processes. 
The numerical labels on the figures denote the thickness of the barrier of the nanostructure and the constant satisfies 
(j = 2e 2 /ha 2 . 



where a is the lattice spacing, e is the electrical charge, fi is Planck's constant, R n is the junction resistance, 
f(u>) = 1/[1 + exp(u/T)] is the Fermi-Dirac distribution, pbulk is the interacting bulk DOS of the barrier 
material, and L = Na is the thickness of the barrier. 
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The resistance-area product versus temperature is plotted in Fig. |4](b). Note the flat regions for low T 
and thin junctions. This is a signal that the transport is dominated by tunneling, but because the gap is so 
small in this system, we do not see an exponential increase in the resistance-area product with the thickness 
of the junction. Instead, it increases with a functional behavior in between that of an exponential increase 
and of a linear increase. We also include a plot of the points where the Thouless energy is equal to the 
temperature. In this case, they are close to, but not exactly at the point where the tunneling behavior crosses 
over to incoherent transport (which has a strong temperature dependence). As U is increased further, this 
separation becomes more readily apparent |3|. We believe this anomalous behavior occurs because the 
system is so close to the critical point of the metal-insulator transition, 

4 Conclusions In this contribution we examined properties of a nanostructure composed of metallic 
leads and a barrier that could be tuned through the metal-insulator transition (described by the Falicov- 
Kimball model). We concentrated on general properties of the thin barrier (where we showed it always 
has a metallic DOS generated by the normal proximity effect with the metallic leads, although the weight 
within this metallic "subband" can be quite small). We also investigated the DOS and transport properties 
of a near critical Mott insulating barrier with U — 5 (the critical value of the transition on a simple cubic 
lattice is (7 C k 4.9 within DMFT). We found the system shows behavior that looks like tunneling, but it also 
has a number of anomalies, the most important being that the resistance does not grow exponentially with 
the junction thickness in the tunneling regime. We also defined a generalization of the Thouless energy that 
reduces to the ballistic and diffusive limits, but can also describe strongly correlated insulators. We found 
the Thouless energy picks up strong temperature dependence in this regime and the point where Eth = 
T determines an approximate crossover from tunneling to incoherent transport, where the resistance is 
proportional to the bulk resistivity of the barrier multiplied by some geometrical factors (and the resistivity 
has strong exponentially activated behavior in T for a correlated insulator). 

There are a number of future directions that are important to consider as one tries to examine models 
closer to experimental systems. First, there will be a charge transfer (electronic charge reconstruction) 1121 
1131 when the chemical potentials do not match, which can have strong effects on correlated systems, 
especially correlated insulators, second, it is useful to consider capacitive effects for these devices, since the 
junction capacitance will play a role in the switching speed, and third, it would be interesting to extend this 
analysis from equilibrium/linear response to nonequilibrium/nonlinear response, where one could directly 
calculate the current-voltage characteristic of the device, and determine the origins of its nonlinearities. 
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